CLNS 01/1763 
CLEO 01-21 



Lifetime Differences, Direct CP Violation and Partial Widths in 
Meson Decays to K^K and tt+tt^ 

(CLEO Collaboration) 
(6 November 2001) 



Abstract 

We describe several measurements using the decays K'^K~ and vr^vr^. 

We find the ratio of partial widths, T (Z)° K+K~) /T vr+vr-), to be 

2.96 lb 0.16 lb 0.15, where the first error is statistical and the second is system- 
atic. We observe no evidence for direct CP violation, obtaining Acp{KK) = 
(0.0 lb 2.2 lb 0.8)% and y4cp(7rvr) = (1.9 ib 3.2 ib 0.8)%. In the limit of no CP 
violation we measure the mixing parameter ycp = —0.012 ib 0.025 ib 0.014 by 
measuring the lifetime difference between — > K'^K~ or vr+vr^ and the CP 
neutral state, — > K~it'^ . We see no evidence for mixing. 
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The structure of the Standard Model has been guided by measurements of mixing and 
CP violation in the neutral K and B meson sectors. The Standard Model predictions for the 
rate of mixing and CP violation in the charm sector are small, with the largest predictions 
in both cases being C(O.Ol), and most predictions being 0(0.001) [Q. Observation of CP 
violation above the 1% level would be strong evidence for physics outside the Standard 
Model. 

The SU(3) flavor symmetry predicts T ^ K+R-) /T (D° ^ Tr+Tr") = 10, while the 
previously measured value is 2.80 ± 0.20 |Q. This deviation is most likely caused by large 
final state interactions. These can also give rise to a large strong phase differences between 
mixing and Cabibbo-suppressed decays that give rise to the same final states [Q. A 
measure of CP violation in these decays, the direct CP violation asymmetry, is proportional 
to the amount of CP violation in the decays and the sine of the strong phase difference. The 
Standard Model suggests that CP violation in these decays is small since the higher-order 
diagrams are suppressed, however, new physics can enhance the rate of CP violation. In 
this paper we present the most precise measurement to date of the ratio of partial widths, 
r [D^ — s> K^K~) /r {D^ — s> TT+TT^) 1^. We also present our search for direct CP violation in 
these decays. 

In the absence of CP violation, the D meson mass eigenstates D12 are also CP eigenstates. 
The decay of a D° to a CP eigenstate, such as K^K~ or tt+tt", has a purely exponential 
lifetime characteristic of the associated mass eigenstate. Therefore, in the limit of no CP 
violation, we can write the time-dependent rate of a decaying to a CP eigenstate, /, as 
R{t) oc exp [— tr ■ (1 — UcpVcp)]^ where CP \ f) = rjcp | /), F is the average D° width, 
ycp = y = Ar/2r, and Ar is the width difference between the two mass eigenstates [||. 
We can measure ycp simply by measuring the ratio of lifetimes of the decaying to a CP 
eigenstate {tcp+) and a CP neutral state such as K'tt^ (r). Then ycp = t/tcp+ — 1. We 
have used r = {tcp+ + tcp-)/2, and assumed that the lifetime difference is small so that 
the K'^n' lifetime distribution can be fit with a single exponential. 

The data were collected using the CLEO II. V upgrade of the CLEO II detector P 
between February 1996 and February 1999 at the Cornell Electron Storage Ring (CESR). The 
data correspond to 9.0 fb~^ of e~^e~ collisions near s/s ~ 10.6 GeV. The detector consisted 
of cylindrical tracking chambers and an electromagnetic calorimeter immersed in a 1.5 Tesla 
axial magnetic field, surrounded by muon chambers. The reconstruction of displaced vertices 
from charm decays was made possible by the addition of a silicon vertex detector (SVX) in 
CLEO II. V. We utilized this improved resolution in previous searches for D^-D^ mixing |Q 
and in measurements of charmed particle lifetimes [|1^. The charged particle trajectories 



were fit using a Kalman filter technique that takes into account energy loss as the particles 



pass through the material of the beam pipe and detector |TT1] . 

The events are selected by searching for the decay chain D*^ -D^tt^, with subsequent 
decays of the to K~^K~, 7r"'"7r~, or K~7r~^. The charge of the slow pion, tt^, from the 
D*~^ decay is a tag of the initial flavor. Additionally, we separate signal from background 
using the energy release in the D*~^ decay, Q = M* — M — M^r, where M* is the candidate 
D*^ invariant mass, M is the candidate invariant mass, and M^r is the pion mass. 

All pairs of oppositely-charged tracks of good quality are used to form candidates 
assuming four particle assignments: K~^K~, K'^n', n~^K~, and vr+vr". The D'^ candidate 
is retained if any of the particle assignments has an invariant mass within 35 MeV of the 
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mass. The D° daughters are constrained to come from a common vertex, and the confidence 
level from this constraint must be greater than 0.01%. A pion candidate with at least two 
SVX hits in both the r-0 and r-z layers is combined with the candidate to form a 
D*^ . The slow pion candidate is refit by constraining it to come from the intersection of the 
beam spot and the projection of the momentum vector. This dramatically reduces the 
mismeasurement of the pion momentum due to multiple scattering in the beam pipe and 
first layer of silicon. The resulting Q distribution has a width of approximately 190 keV. 
The candidate is retained if the confidence level for the refit is greater than 0.01%, Q is 
less than 25 MeV, and the D*^ momentum is greater than 2.2 GeV/c. Finally, we require 
I cos^* I < 0.8, where cos^* is the angle in the rest frame between a daughter and the 

direction in the lab frame. The signal is flat in cos^^*, while the backgrounds are highly 
peaked at | cos^^* \^ 1. Particle identification using specific ionization is not required since 
the different mass hypotheses are separated by greater than 8.5 standard deviations. 

The partial width measurements are obtained from binned maximum likelihood fits to 
the Q distribution of the D*^ decay. We fit in bins of momentum to eliminate potential bias 
due to mismodeling of the D*^ momentum spectrum in Monte Carlo. The finite statistics 
of the fitting shapes are included in the statistical uncertainty of the fit. The signal shape is 
taken from the K'-r^ data sample, while the background shape is determined from Monte 
Carlo. All of the modes have approximately the same signal shape since the Q resolution is 
dominated by multiple scattering of the slow pion. We also fit the mass distribution as 
a check. 

We first fit the Kti data outside of the signal region to obtain the background nor- 
malization. To obtain i^jrvr = T {D^ tt^tt^) /V {D^ K~7r~^) we fit the Q distributions 
for the ratio of signal yields between the tttt and Ktt channels, and for the normaliza- 
tion of the background, where we have used the signal shape and background parame- 
ters determined from the Kn data and Monte Carlo samples, respectively. To obtain 
Rj^j^ = T{D° ^ K+K-) /r (1}° K-n+) we fit the Q distributions as we did for R^^, 
however we add an additional component from pseudoscalar-vector decay (PV) background, 
where the shape is taken from Monte Carlo and the normalization is allowed to fioat. The 
PV background is primarily from D° — > K~ , p+ — > 7r''"7r° where the 7r° is nearly at rest. 
This background forms a broad peak in Q. The PV background is negligible in the tttt and 
Ktt final states. 

In order to maintain statistical independence, we use two different sets of Monte Carlo 
events. One sample is only used to determine the fitting shapes. We fit the data and the 
second Monte Carlo sample simultaneously to correct for small differences in acceptance 
between the normahzation and signal modes. The results of the fits are Rrk — 0.1037 ± 
0.0038 and R^^ = 0.0355 ± 0.0017 from approximately 20,000 X-7r+, 1900 K+K', and 710 
7r"'"7r~ events. 

The systematic uncertainty due to the fitting shapes is assessed by performing a series 
of fits using different assumptions for the background and also several fits to the mass 
distribution. We estimate systematic uncertainties of 0.0017 and 0.001 due to the fitting 
shapes in the KK and vrvr modes, respectively. We vary the bin sizes, Q fit range, Q signal 
region and candidate mass requirement, and form a combined systematic uncertainty of 
0.0005 due to these variations. 

We also estimate systematic uncertainties associated with some of the event selection 



5 



requirements by doing the analysis without those requirements. The variations we observe 
are 0.0009 in Rkk and 0.00095 in i?7r7r from removing the vertex confidence level requirement, 
and 0.00032 in Rrk and 0.00016 in Rt^t^ from removing the track quality requirement. 

We use the Kt: data sample to study the effect of any mismodeling in the simula- 
tion of the fragmentation and the detector acceptance. For the fragmentation model- 
ing we estimate a systematic uncertainty of 0.0014 for Rkk and 0.0005 for i?^^. We 
obtain relative corrections and uncertainties due to mismodeling of the detector accep- 
tance of (-2.4 ± 1.1)% for Rkk and (+2.4 ± 2.7)% for R^^. We apply these correc- 
tions and sum all of the systematic uncertainties in quadrature to obtain the final re- 
sults Rkk = T{D° ^ K+K') /T K-n+) = 0.1040 ± 0.0033 ± 0.0027 and R^^ = 
r {D^ TT+TT-) /r (D^ K-n+) = 0.0351 ± 0.0016 ± 0.0017, where the first error is sta- 
tistical and the second is systematic. These results are the most precise determinations of 
Rkk and R.„n to date [Q. 

We can combine the results, accounting for cancellations and correlations among the 
uncertainties to calculate Rkk/Rtt-k = 2.96 ± 0.16(stat) ±0.15(syst). This result agrees with 
the world average value of 2.80 ± 0.20 gj. 

We can use the same procedure to search for the direct CP asymmetries 

where / can be K^K~ or 7r"^7r~. The charge of the slow pion from the D*^ decay serves as 
an unbiased tag of the fiavor since charm quarks are produced in quark-antiquark pairs 
at CESR and fragmentation and the D* decay are strong processes, which conserve CP. 

We measure the CP asymmetry in the same manner as the partial width analysis de- 
scribed above apart from the following changes. The K^K~ and tt+tt^ data are separated 
into and samples based on the charge of the slow pion. However, we still normalize by 
the entire Ktt sample to eliminate possible bias from any asymmetry in — * K~ti^ decay. 
The D*^ momentum requirement is altered to be greater than 2.0 GeV/c since acceptance 
differences between modes are no longer an issue. The candidate mass requirement is 
tightened to ±15 MeV of the nominal mass, which reduces the backgrounds by about a 
factor of two. 

We fit the data in the same manner as in the partial width analysis, modified as described 
above. The KK and tttt Q distributions and fit results are shown in Figure From the 
fits we find 1512 ± 47 ^ K+K_- events, 1511 ± 47 ^ K+R- events, 579 ± 26 
TT+TT" events, and 557 ± 26 ^ tt+tt" events, and obtain A^^ = 0.001 ± 0.022 and 
A^^'p = 0.020 ± 0.032. 

The sources of possible systematic error for the CP asymmetry measurement are the 
shapes used for fitting and a charge dependent slow pion acceptance. To assess the systematic 
uncertainty from the fitting shapes we perform fits in which we vary the candidate mass 
window, remove the vertex confidence level requirement, vary the width of the i^Tr signal 
region and the Q fit region, alter the number of bins, and split the Kir sample into two 
according to the charge of the associated slow pion and fit the two samples separately. We 
use 1/2 of the largest variation in each case, and then sum them in quadrature to obtain a 
systematic uncertainty due to the fitting shape of 0.0068 for A^p and 0.0069 for A^p. 
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FIG. 1. The D* Dtts Q distributions for a) D° K+R- and D° K+K~ candidates 
and b) — > vr+vr^ and vr^vr^ candidates. The points are the data and the histograms are 

the background fits. 



A difference in slow pion acceptance for positive and negative pions can come from a 
number of different sources. The interaction cross section of pions with matter is different 
for positive and negative pions. We use the known composition of the CLEO detector and the 
interaction cross sections to calculate the induced asymmetry as a function of momentum. 
We find that the bias to the asymmetry is less than 0.2%. We use the pions from decays to 
search for a momentum-dependent charge bias in pion acceptance. We select the pions from 
Kg decay similarly to the method used to select the slow pions from D*~^ decay. We compare 
the observed difference between the momentum spectrum for the positive and negative legs 
of the K^, over the region of slow pion momenta from D*~^ — >■ D^vr"*" decay, to estimate the 
acceptance difference for positive and negative pions to be less than 0.07%. 

We have looked for a momentum-independent charge bias in track finding by generating 
single track Monte Carlo randomly distributed in 6, (p, and momentum, between and 3 
GeV/c. We see no significant bias, and limit the momentum independent acceptance bias 
to be less than 0.48%. We translate these limits on acceptance differences and track finding 
biases into limits on our observed asymmetry based on the statistics of our observed data 
sample. 

Charm quarks are expected to be produced with a small forward-backward asymmetry 
in e~^e~ annihilations at ^/s ~ 10.6 GeV due to the interference between the photon and 
The center of the luminous region was not exactly at the center of the detector, so 
this, coupled with the forward-backward asymmetry, induces an acceptance asymmetry. 
From a study of the K~^tt~ data and Monte Carlo samples we find an acceptance bias of 
0.014 ± 0.014%. We correct for the bias and assign the statistical error as a systematic 
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uncertainty. 

Summing all of the systematic uncertainties in quadrature and applying the correction 
mentioned above we arrive at the final result of A^j^ = (0.0 ± 2.2 ± 0.8)% and A^^'p = 
(1.9 ± 3.2 ± 0.8)%. We see no evidence of direct CP violation in these decays. This is the 
most precise measurement of these CP asymmetries to date |Q,|T2 . 



As noted earlier we can measure the normalized mixing parameter ycp by measuring 
the lifetime ratio between K^tt~^ and decay to a CP eigenstate, such as K^K~ 

or TT+TT": ycp = t/tcp+ — 1. In the limit of no CP violation in the D meson sector ycp 
is equivalent to y. We use the same data sample described above, using the decay length 
and momentum to determine the proper decay time. We modify the event selection criteria 
slightly for this analysis. We require the candidate momentum to be greater than 2.3 
GeV/c. We tighten the requirement on the vertex confidence level of the D° candidate to 
be greater than 0.1%. Furthermore, we place an extra requirement on the data: the 
candidate masses obtained with the three other particle assignments to the two daughters 
must be more than four standard deviations away from the nominal D° mass. 

We select events with a Q value within 1 MeV of the nominal value and fit their candidate 

mass spectrum with a binned maximum likelihood fit to the sum of two Gaussians for the 
signal, constrained to the same central value, and a first order polynomial for the background. 
The data and fit results are shown in Figure ^ The fit values are converted into a mass 
dependent probabihty for signal and background and are used as an input to the lifetime fits. 
The other inputs to the lifetime fits are the measured proper decay time and its calculated 
uncertainty. For the KK and tttt samples we fix the ratio of areas and the ratio of widths of 
the two Gaussians to the values determined in the Ktt fit. We perform the fits for candidate 

mass over the range 1.825 to 1.905 GeV, and use all of these events in the lifetime fits 
described below. 

For the signal portion of the probability distribution function for the lifetime fits we 
constrain the candidate mass to a fixed value, which gives us a better measurement of 
ycp- The value we constrain to is the weighted average of the mass determined from 
the Kit, KK, and tttt events, where each is corrected by an offset determined from Monte 
Carlo. This offset is simply the difference between the input and measured mass for each 
channel in the Monte Carlo. The offsets are +0.15 ± 0.02 MeV (Kn), +0.27 ± 0.05 MeV 
[KK), and +0.10 + 0.09 MeV (tttt). The mass constraint introduces a systematic bias in the 
lifetime measurement, which cancels for ycp which only depends on the ratio of lifetimes. 

The candidate proper decay time, t, is given by 

. V dec ^prodJ ' P 
t = m ■ , 

I P I 

where r^ec and p are from the candidate vertex fit. We determine rprod using e+e" 
qq [q = udsch) events from sets of data with integrated luminosity of several pb~^. The 
extent of the luminous region has a Gaussian width of approximately 10 /xm vertically, 300 
/im horizontally, and 1 cm along the beam direction [0. The resolution on the decay 



point is typically 40 /xm in each dimension. The resolution in t is typically at = 0.4 in units 
of lifetimes. We determine the proper decay time in the three dimensions separately, and 
combine them to arrive at the best estimate of t and ctj. 
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FIG. 2. The mass distribution for K^K' (left) and vr+TT (right) candidates. 

The curves are the results of the fit discussed in the text. 
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TABLE I. Summary of the lifetime fits. The parameters are those described in the text, where 
/mis is the fraction of signal in the second and third Gaussian contributions and fTmis is the width of 
the second Gaussian. Note that we have constrained the candidates to a mass of 1.86514 GeV, 
the Monte Carlo corrected weighted average of the KK, vrvr, and Ktt data. This mass constraint 
introduces a systematic bias in the lifetime measurement, which cancels for ycp which only depends 
on the ratio of lifetimes. This technique yields the smallest uncertainty in ycp, but is not optimal 
for measuring the absolute lifetime. 



Parameter 


Ktt 


KK 


vrvr 


Number of signal 


20272 ± 178 


2463 ± 65 


930 ± 37 


Tsig (ps) 


0.4046 ± 0.0036 


0.411 ±0.012 


0.401 ±0.017 


Background frac. (%) 


8.8 ±0.2 


50.7 ±0.7 


29.1 ± 1.3 


Background life frac. (%) 


81.0 ±4.8 


85.7 ±2.9 


32.2 ±7.5 


7"back (ps) 


0.376 ± 0.030 


0.436 ± 0.020 


0.56 ±0.15 


/mis (%) 


3.8 ±0.9 


Fixed 


Fixed 


0-mis (ps) 


0.590 ± 0.079 


Fixed 


Fixed 



We fit the lifetime distribution using an unbinned likelihood method. The signal prob- 
ability distribution function (PDF) consists of an exponential convolved with a resolution 
function, composed of the sum of three parts, based on a simple, yet robust, physical model. 
For most events the calculated covariance matrix for the daughters is assumed to be 
correct to within a global scale factor, with a Gaussian resolution function of width S ■ at- 
The scale factor, S, accounts for any common mistake in the covariance matrices, as would 
be present from a deficiency in the detector material description. A few percent of the events 
have one or more particles that have undergone a hard scatter, rendering the extrapolated 
vertex errors virtually meaningless. We model the contribution from these events with a 
single Gaussian whose normalization and width are allowed to float in the fit. For a very 
small fraction of events the vertex location is extremely mismeasured. These events have 
a nearly flat distribution in lifetime. We model this contribution with a broad Gaussian, 
assigning a fixed width of 8 ps. The normalization of this contribution is allowed to float 
in the fit. The signal PDF is multiplied, on an event-by-event basis, by the mass-dependent 
signal probability from the candidate mass fit. 

The background lifetime distribution contains two pieces: a prompt piece and a piece with 
non-zero lifetime. The component with non-zero lifetime comes from partially reconstructed 
charm decays. We model this component with a single exponential where the lifetime is 
another parameter of the fit. We expect the fitted value of the background lifetime to be 
consistent with the lifetime. The relative amount of background with and without lifetime 
is also allowed to float in the fit. Both sorts of background are convolved with a resolution 
function that is modeled in the same manner as the signal, but with an independent set of 
parameters. The background PDF is multiplied by one minus the signal probability from 
the mass fit. 

The fit results for all events included in the fit are shown in Figure ^ and given in Table |. 
Small corrections to the lifetimes are computed by comparing the generated and measured 
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KK Proper Time (ps) njt Proper Time (ps) 

FIG. 3. The proper time distribution for all — > K^K~ (left) and — > tt+tt" (right) 
candidates included in the fit. The curves are the fit results discussed in the text. 
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values in a Monte Carlo analysis on a fully simulated sample, including backgrounds, cor- 
responding to roughly ten times the data sample. These corrections are 0.0006 ± 0.0040 ps 
in K+K-, -0.0011 ± 0.0015 ps in K-ti+, and 0.001 ± 0.0058 ps in vr+vr-. Applying these 
corrections we obtain y^j^ = -0.019 ± 0.030 ± 0.010, = 0.005 ± 0.046 ± 0.014, and 
combining them in a weighted average we calculate ycp = —0.012 ± 0.025 ± 0.009, where 
the second error is from the Monte Carlo statistics. 

We check the data for bias in several different parameters. We plot the fitted value of 
ycp versus azimuthal angle, polar angle, date the data were collected, momentum of the 
candidate D^, cos 6'*, and confidence level of the vertex constraint. We find no significant 
biases in any of these distributions. 

The kinematics of Kir, KK, and tttt decays are slightly different due to the differ- 
ent amount of kinetic energy released. This will result in the signal resolution functions 
being slightly different. We have constrained all of the signal resolution functions to be the 
same. Studying this effect in Monte Carlo and data we estimate the following systematic 
uncertainties: 0.007 for KK, 0.003 for tttt, and 0.005 for the average. 

We study the effects of background shape mismodeling by varying the amount and compo- 
sition of the background. We perform these in data and Monte Carlo and estimate systematic 
uncertainties of 0.008 for KK, 0.011 for tttt, and 0.008 for the average. 

We study the effect of our treatment of the proper time outlier events, which we have 
modeled with a wide Gaussian of fixed width. We vary the value of the width used in 
the wide Gaussian and also eliminate the wide Gaussian from the resolution function and 
impose a maximum proper time limit instead. From these studies we estimate systematic 
uncertainties of 0.002 for KK, 0.001 for tttt, and 0.002 for the average. 

We investigate the bias introduced by constraining all of the events to the same 
mass by removing this constraint. We take the difference between the constrained and 
unconstrained fits as a systematic uncertainty: 0.005 in KK, 0.005 in vrvr, and 0.005 in 



the average. Length scale uncertainties have been studied previously by CLEO [|10[ and 
contribute negligible uncertainty to ycp- 

Summing all of the listed systematic uncertainties in quadrature, including the Monte 
Carlo statistics, we obtain the final results |/^^ = -0.019 ± 0.029 ± 0.016, = 0.005 ± 
0.043 ± 0.018 and combining the two results we obtain ycp = —0.011 ± 0.025 ± 0.014, which 
is consistent with zero and with the average of previous measurements [III], ycp = 2.9 ± 1.4. 

In summary, we have used the CLEO II. V data set to obtain the world's most precise 
measurements of Rkk = T {D° ^ K+K-) /T K-'k+) = (10.40 ± 0.33 ± 0.27)% and 

R^^ = r (DO ^ TT+TT-) /r {D^ K-n+) = (3.51 ± 0.16 ± 0.17)%, and the direct CP asym- 
metries A^l^ = (0.0 ± 2.2 ± 0.8)% and A^'^ = (1.9 ± 3.2 ± 0.8)%. We have also performed a 
competitive measurement of the normalized mixing parameter ycp = —0.012 ±0.025 ±0.014. 
In all cases the first error is statistical and the second is systematic. Our partial width mea- 
surements are consistent with the previous world average, we see no evidence for direct CP 
violation in Cabibbo-suppressed decays, and we measure a value of the mixing parameter 
ycp consistent with zero. 
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